Generation of macroscopic quantum-superposition states by linear coupling to a bath 
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We demonstrate through an exactly solvable model that collective coupling to any thermal bath induces 
effectively nonlinear couplings in a quantum many-body (multi-spin) system. The resulting evolution can 
drive an uncorrected large-spin system with high probability into a macroscopic quantum-superposition 
state. We discuss possible experimental realizations. 
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Introduction: The strive for ultrafast and ultrapowerful 
data processing by quantum information techniques com- 
monly relies on quantum entanglement (QE) that is in- 
duced by (direct or indirect) interactions among the par- 
ticles in the system Yet this QE is extremely vulner- 
able to the environment: its fragility exponentially mounts 
with the strength of the system-environment (bath) cou- 
pling and the number of particles Alternatively, QE 
is realizable by collective (multipartite) dissipation, i.e., 
real-quanta exchange via the bath within the rotating wave 
approximation |16fl. In both forms of the QE, the bath ef- 
fects are commonly treated within the Markov approxima- 
tions iHIBlS]- Here we take an alternative approach and 
show through an exactly solvable model that QE can arise 
spontaneously from dispersive, nondissipative interactions 
(virtual-quanta exchange) among particles via the bath, a 
process unexpectedly revealed by going beyond the fore- 
going standard approximations. 

Such bath-mediated dispersive interactions result in ef- 
fectively nonlinear couplings, although we assume the 
system-bath coupling to be linear, as usual They are 
shown to grow up to an asymptotic value that is reached 
beyond the non-Markov (memory) time of the bath. The 
resulting unitary evolution can allow with high probabil- 
ity for the formation of entangled, macroscopic quantum 
superposition (MQS) states in atom or spin ensembles col- 
lectively coupled to any bosonic bath, at zero (vacuum) or 
nonzero temperature. The number of particles in the MQS 
depends on the spectral response of the bath and its temper- 
ature. This result generalizes the notion of QE via disper- 
sive single -mode interactions in ion-traps |l8|| to arbitrary 
bath spectra and diverse scenarios. 

Collective dynamics of ensembles of atoms and spins are 
among the few well-studied manifestations of quantum be- 
havior on macroscopic scales. This comes about because 
their quantized collective dynamics can be mapped onto 
that of an object in an eigenstate of angular momentum 
(spin) L with large eigenvalues. Behavior of this kind is 
exhibited by spin-polarized ensembles in solids Q or by 
atomic ensembles with large pseudospin that collectively 
emit and absorb photons lUOLIllll . Alternatively, large spin 
characterizes macroscopic atomic ensembles that are en- 



tangled via interaction with a common light source 111211 . 

Model and dynamics: We consider an ensemble of 
non-interacting spins or atomic two-level systems (TLS) 
that are identically, linearly coupled to a bosonic (oscilla- 
tor) bath. The model is described in the collective basis by 
the many-body Hamiltonian 

H = Hs + Hb + Hj, 
Hs = i-^tLx 

HB = Ek^kb%, Hj = L^Z,m{h + bl). (1) 

Here bl , bk are the creation and annihilation bosonic op- 
erators of the /c-th bath mode, and t]^ the corresponding 
(k — mode) coupling rates. The collective spin operators 
in Hs and Hj are Li = '^y.(T].{i = x,y,z), the com- 
ponents of the total spin L of the ensemble, al being the 
TLS operators (Pauli matrices) for the k-th spin-1/2 parti- 
cle. The rotating wave approximation has not been made 
here. Since H commutes with = J2i ^1' the bath inter- 
acts separately with each subspace of the system labelled 
by the total-spin value /. The I values range from to y. 



if is even, and 



to ^, if TV is odd. 



Hence, it is suffi- 
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cient to study the interaction of the bath with a multi-level, 
(2/ + l)-dimensional system. 

In general, the dynamics generated by Eq. ([T]l is insolv- 
able, because of the non-commutativity of and L^. In 
order to circumvent this difficulty we prepare the system 
in an eigenstate of L^, = J2k^k superposition of 
eigenstates) and then switch off u^. This removes Hs in 
Eq. ([ll. The Hamiltonian dynamics then becomes exactly 
solvable for any bosonic bath. For a given spin-sector / of 
the ensemble, one can then write a closed-form equation 
for the time-evolution operator of the system and the bath 
(see lll3n for more details), given by 



ii/(t)L^ + L,5](a,(t)6l-a*(t)6, 



Ui{t) = exp 



where the coupling spectrum of the bath (see below) deter- 
mines the functions 



(2) 



f{t) = j^r]l{u}kt - sinujkt)/ujl, ak{t) = r]k 
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We thus obtain a striking exact result: the bath-induced 
evolution is driven by both (linear) and {nonlinear) 
terms. The linear terms cause decoherence as expected. 
The new term f(t)Ll in Eq. ^ occurs for multipartite 
systems where = X + J2i=;ij '^i^j' absent in the 
single particle case where = = I. As discussed 
below, it may be interpreted as a collective analog of the 
Lamb shift Ill4h or frequency pulling of each spin by all 
others via the bath (virtual quanta exchange). 

We focus on a class of initial states of the spin ensemble 
that can be expressed as a weighted sum of density matrices 
with values of /, p(0) = J2i ^iPii^)- Under the above 
dynamics each component pi evolves separately. The state 
of the system (initially uncorrected with the bath), at any 
later time, is found upon tracing over the bath degrees of 
freedom: 




Pit) = E 



itf(t)Ll 



tr(t)(rn-m')^ 



\m){m'\ 
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(4) 

where T{t) = j J2k vli^ ~ cosujkt)/Lo1 is the decoher- 
ence rate (discussed below) and m, m' , ranging from +/ to 
-I, label the states ifToll . 

In keeping with the procedure leading to Eq. (|2|, at time 
t = each spin is prepared in a superposition of its en- 
ergy (o"|) eigen states The total system is then initially in a 
product of such superposition states which is uncorrected 
(unentangled) among the individual spins. The initial un- 
corrected state of the spin ensemble (all spins identical), 
can then be written as p(0) = |^(O))('0(O)|, where 



IV^(o)) 



.4,) = |V'i><»IV'2>---|V'jv>, IV'fc 



COS -I t)+sin-e" 
(5) 



6 and (f) are the usual Bloch-sphere angles. This state is an 
eigenstate of the collective spin operator L ■ h and n is the 
unit vector corresponding to the angles and (f). 

The off-diagonal terms in Eq. dU decay exponentially at 
the rate r(t)(77T, — m'Y, rendering multipartite coherence 
(entanglement) vulnerable to decoherence. Let us, how- 
ever, first consider that T{t) is negligible. Then, under the 
nonlinear term f{t)Ll in Eq. Q, the initial uncorrelated 
state (|5]l evolves at prescribed times (see below) into an en- 
tangled macroscopic quantum superposition (MQS), state 
(analogoulsy to evolution under the nonlinear Kerr Hamil- 
tonian H): 



Figure 1 : (Color online) Bath-induced formation of a macro- 
scopic quantum-superposition (MQS) state, (a) A collective 
(large) spin naturally evolves in the bath into a MQS state. The 
system is a spin ensemble composed of A*' = 50 particles cou- 
pled to an Ohmic bath. The cut-off frequency ujc is chosen 
such that the collective coupling of the spin ensemble to the 

bath, rj = \J^k ~ O.OOScjc. (b) The coupling spectrum 
of an Ohmic bath, (c) The time-dependent functions responsi- 
ble for the nonlinear Lamb-shift (/(f)) and decoherence (r(t)) 
dynamics of the system, (d)-(e) The absolute value of the den- 
sity matrix elements, \pmm'{t) \ in the Lx basis (Eq. (O) at var- 
ious times. The initial state of the ensemble is a spin coher- 
;n|t^state, \d = 7r/4, = 0), in Lx basis. The presence of the 
±N/2) {^^N/2\ off-diagonal elements (left and right most peaks) 
signifies the formation of such a state with high fidelity. 



negligible decoherence (r(t) ~ 0), into a MQS which is 
simultaneously oriented along the +x and —x directions 
with vr /2 relative phase 
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^/^|^,(/.)+e^'^/^|^-^, (/))]. (6) 



If, for example, = |f ,0), then Eq. (|6]l corre- 

sponds to a superposition of the state with / = N/2, m = 
N/2 in the Lx basis, where all the spins are oriented along 
+x direction, and the state with / = N/2,m = —N/2 
in Lx basis, where all the spins are oriented along —x di- 
rection. The interaction with the bath in Eq. (|4]l thus trans- 
forms the initially uncorrelated state, | = | f , 0), under 



This is a macroscopic GHZ-like state in which all the N 
spins are maximally entangled. As is known, the parity 
of N plays a role 10, IH: only even N strictly result in 
GHZ states, although for large N the deviation from per- 
fect GHZ states for odd N is negligible. 

The MQS will form according to Eq. ([3]), dU at times 
when tf{t) = {2n -\- l)'/r/2 (n being an integer) The 
earliest formation time of the MQS is then 



Tmqs = t 



vr 



2/(t) 



(8) 



The time at which such a state forms is independent of A^. 
Other entangled states ("Schroedinger kittens") will form 
at intermediate times iB flSll . 
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Since the bath spectrum is continuous, the decoherence 
r(t) 7^ is unavoidable and the formation of perfect MQS 
states is generally not possible. Yet, the formation of high- 
purity states i.e., the accomplishment of Eq. (|6]l at Tmqs 
(Eq. dS])) with high probability may still be possible, as ar- 
gued below. The condition for the quantum-superposition 
to survive decoherence is, from Eqs. (|4]l-(l8]l, that the MQS 
forms faster than it decays: 



Tmqs^N^ < 1- 



(9) 



where TN^ is the upper limit of the time-averaged decay 
rate T{t){m — m'Y of the off-diagonal (coherence) ele- 
ments in Eq. Q. 

Dependence of MQS formation on bath spectrum and 
temperature: What determines whether the condition Eq. 
(|9]l for MQS formation can be satisfied ? This requires that 
f[t) strongly exceeds T{t) at the time Tmqs- The func- 
tions f{t) and r(t) in Eqs. (|3]l, Q, dSJ are then expressed 



as 
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cot — sin ujt 



^ , , 1 — cos Ujt 
Gt{uj) , (10) 



where Gq{uj) and Gt{^^) are, respectively the zero- 
temperature and finite-temperature coupling spectra of the 
bath. The temperature -dependent coupling spectrum of the 
bath is defined as fl [3 Gt{lo) = J2k^ll duj^i^ - 
u}k)cothf3u}, where f3 is the inverse temperature of the 
bath. It is seen from Eq. (fTOll th at f{t) is related 
to the bath-induced Lamb shift UM (the real part of 
the zero-temperature bath susceptibility) and T(t) to the 
bath-induced decoherence rate (the ima gina ry part of the 
temperature-dependent susceptibility |[3l Il4ll ). Hence, the 
feasibility of condition (|9]l, which requires f{TMQs) ^ 
T{tmqs)' is determined by the bath coupling spectrum and 
temperature. 

It is advantageous for the satisfaction of condition (|9]l 
that the Tmqs exceeds the non-Markov time scale for the 
following reason. The decoherence rate T{t) is time- 
dependent in the non-Markov regime of t < tc, where 
tc, the correlation (memory) time of the bath is the in- 
verse width of Gt{^^) At sufficiently low temper- 
atures, r(t) is drastically reduced in the Markovian limit 
(t 3> tc) as opposed to its fast initial non-Markovian in- 
crease: r{t < tc) > T{t oo) = Fm (Fig. 1(c), 
2(a)). This comes about since T{t) that initially has con- 
tributions from all the bath modes, J GT{uj)du}, subse- 
quently decreases, as the bath mode-oscillators go out of 
phase in the Markov regime lll6ll . By contrast, f{t) in- 
creases in the course of the transition from the non-Markov 
to the Markov regime, where it settles at its long-time 
value \f{t oo)| = |/a/| » \f{t < tj]. Since the 
MQS-state formation time Tmqs is typically longer than 
the bath correlation time tc, MQS encounters a much lower 



r ~ T]\j, and much higher / ~ f^, than those encoun- 
tered in the non-Markov regime, thereby helping satisfy 

m- 

At sufficiently high temperatures we attain the regime 
Fm > /m where the formation of MQS is not possible. 
For an Ohmic bath the condition /m > Fm is satisfied 
when the cutoff frequency (energy) is larger than the 
thermal energy, u!c > ksT [H = 1) lfl3ll . 

These trends can be seen from Fig. 1(c) for an 



Ohmic bath with coupling spectrum Gq = aue 



V-^c and 



Fig. 2(a), for a Lorentzian bath with coupling spectrum 

2 

Go{u!) = a ^2^("l^„i2 ■ Different bath spectra having 



the same width, i.e., the same inverse correlation time 1/tc, 
may still have different fj^i values and different Fm in the 
long-time Markov limit and hence yield MQS with differ- 
ent purities at tmqs according to (Eq. ©). 

Spin ensemble coupled to a phonon bath: As an exper- 
imentally feasible scenario, consider a localized aggregate 
of N weakly interacting spin- 1/2 particles (e.g. electron 
spin ensembles in fullerenes ||9l[l3l or quantum dots ifTsIl ) 
that undergo dephasing via identical coupling to phonons 
in a lattice, of longer wavelength than the aggregate size 
. Hence, they conform to the model of Eq. ([T]). For an 
Ohmic phonon spectrum with the Debye cutoff loj^ we find 
|[l3ll /m = olujd, and Ff,[ = a//?, where j3 is the inverse 
temperature and a < 1 is a dimensionless coupling con- 
stant. Condition (|9]l may be satisfied at temperatures below 
ImK, for /m ~ u;_d ~ lO^i^z » Fm > MHz to ob- 
tain a quantum-superposition (GHZ state) with A'^ ~ 300 
spins with high fidelity (Fig. 1) at Tmqs- 

Atomic ensemble coupled to a cavity: Another exam- 
ple we consider is a non-interacting atomic gas coupled to 
a single-mode cavity lfl9ll (photonic bath). Whereas for 
isotropic spin ensembles the quantization axes in (1) are 
arbitrary, this is not the case for two-level atoms (TLS) 
i.e., pseudospins lfl0l - [l2ll . There the energy splitting of 
Hs levels is represented by Hs = ^j^zLz- Correspond- 
ingly, L^— coupling of atoms to a photonic bath causes 
pure cooperative dephasing without population exchange, 
while Lj; coupling to the bath causes cooperative popu- 
lation exchange or relaxation. In the latter case, we have 
to initialize the system in a L^-eigenstate and then bring 
each TLS to degeneracy i.e. set oo^ = by applying Zee- 
man shifts. Once the (two-level) atoms are prepared in 
degenerate Zeeman states {Hs = 0) one can induce 
coupling between the cavity and the atomic ensemble by a 
two-photon Raman process . The collective nonlinear evo- 
lution can thus generate MQS of the atomic ensemble. For 
finite cavity linewidth, the cavity acts as a Lorentzian bath, 
which results in dynamics similar to Fig. 2(a). The deco- 
herence rate T^v/ is determined by the zero-frequency cou- 
pling strength to the cavity Gt(0), or may be induced by 
an external low-frequency noise. This may limit the size of 
MQS allowed by Eq. Q to ~ 100, for an atom-cavity 
coupling strength r] ~ IMHz. 
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Figure 2: (Color online) Interaction with a cavity mode. As 

in Fig. 1, for a cavity with Lorentzian coupling spectrum with 
width bJc centered at lOojc. (a) The f{t) and V{t) functions, (b) 
The formation of macroscopic quantum-superposition for cavity 
with Lorentzian lineshape at tmqs = 10^/a;c. 

Conclusions: We have used an exactly solvable model 
to reveal the unexpected entangling dynamics of a system 
with large angular momentum that is linearly, collectively 
coupled to a thermal bosonic bath. The intriguing con- 
sequence is that a commonly occurring finite-temperature 
environment may naturally induce rather than impede the 
formation of macroscopic quantum-superposition (MQS) 
states. Such counterintuitive bath-induced effects change 
our fundamental perspective of non-classicality in open 
quantum systems by identifying a broad class of natural 
entanglement. On the applied side, the feasibility of high- 
fidelity entangled states with N > 100 may be a starting 
poin t to the advancement of one-way quantum computing 

Equations ©, 0, ©, (|9| are the main results of this 
exact analysis. They imply novel dynamic features in open 
multispin systems at low temperatures: (i) Non-linear (L^) 
terms induced by a bosonic bath can dominate the evolu- 
tion of multi-spin systems. By contrast, for a single spin- 
1/2 (two-level) system, -^^^(j,) = '^l(x) ^ 
linear term only yields an overall phase and does not af- 
fect the dynamics, (ii) While the non-Markovian bath dy- 
namics affects the squeezing of an initial coherent state, 
its Markovian dynamics governs the formation of a high 
purity MQS. (iii) The design of a bath-coupling spectrum 
that determines a high ratio of the Lamb shift fm to the 
decoherence rate Tm, can play a crucial role in allowing 
the formation of such macroscopic quantum- superposition 
states. (iv)Their purity may be further improved using 
quantum control techniques to dynamically modulate the 
system levels 1I21I1 and thus modify the bath effects, (v) 
Another intriguing consequence of the analysis is that only 
a bosonic bath induces such nonlinear dynamics, whereas 
a fermionic bath generates also higher-power nonlineari- 



ties which could be unfavorable for the generation of high- 
fidelity cat states. 
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